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Abstract 

The non-hnearity of Einstein's equations makes it possible for small-scale 
matter inhomogeneities to affect the Universe at cosmological distances. We 
study the size of such effects using a simple heuristic model that captures the 
most important backreaction effect due to nonrelativistc matter, as well as 
several exact solutions describing inhomogeneous and anisotropic expanding 
universes. We find that the effects are 0{H'^P /c^) or smaller, where H is the 
Hubble parameter and I the typical size scale of inhomogeneities. For virialized 
structures this is of order v'^/c^, where v is the characteristic peculiar velocity. 



1 Introduction 



The Friedmann-Robertson- Walker (FRW) homogeneous solution to Einstein's equa- 
tions gives an excellent approximation to the large-scale structure of our universe (for 
a survey, see e.g. PP). This seemingly uncontroversial assumption is less obvious than 
it appears because the distribution of matter-visible or otherwise-is inhomogeneous 
already at ~ 10 Mpc scales. At that scale the density contrast 5p/ p becomes 0(1) 
and the non-linearities of Einstein's equations makes it conceivable that this may 
affect some of the properties of the universe even at cosmological scales. 

This possibility has been examined in many papers. Extreme effects were advo- 
cated in [2j, which argued that primordial CMB inhomogeneities coupled with the 
non-linear nature of Einstein's equations could account for the late-time acceleration 
of the universe without the need for dark energy. Convincing arguments against this 
were advanced in j^l E] , yet the less extreme claim that inhomogeneity can affect the 
very large-scale behavior of the universe is worth investigating. Indeed, even if these 
effects turn out to be small, they may still be relevant for interpreting data of the 
next generation of cosmological probes. For instance the study of Ref. [Sj suggests 
that the backreaction effect of inhomogeneities at the Hubble scale H is of the order 
of 10"^ 

When studying the effect of inhomogeneities, one must be careful in not mistaking 
the onset of large matter inhomogeneities with the breakdown of the linear approxi- 
mation for the gravitational field itself. Intuitively, the difference is that gravity is so 
weak that a linear equation for the gravitational field holds everywhere outside black 
holes, so, in particular even when 5p/ p 3> 1 (as on the surface of the Earth). Yet 
making this intuition into a precise statement is difficult. One would need a system- 
atic expansion that takes into account matter non-linearities to all orders, but treats 
the gravitational backreaction perturbatively. That is, one would need an appropriate 
perturbative series in G5p ^. 

We will study this problem heuristically in section 2, by presenting a "paradox" 
that arises already at second order in the G5p expansion, together with its resolution. 
Namely, we shall consider point-like particles of mass m, distributed on a regular 
lattice of size /. At scales larger than /, one would expect to find a uniform FRW 
solution with Hubble constant = SvrGm/S/'^. Yet as we will see the backreaction 
of the point-like sources is formally infinite. It becomes finite and of order PH^ only 
after an appropriate renormalization of the mass m. Physically, this renormalization 
arises because one cannot separate a "bare" mass from its gravitational energy. 

Section 2 is heuristic in that the only backreaction term kept there is the Newto- 

is the Newton constant and everywhere in this paper c — 1. 
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nian gravitational energy. This is not a rigorous procedure, because effects of similar 
size are ignored. To better study backreaction effects we proceed in sections 3 and 
4 to study simple exactly soluble cosmologies: an array of equally spaced parallel 
two-dimensional dust walls in section 3, and in section 4 solutions for arbitrary ar- 
rangements of parallel cosmic strings and the "Swiss Cheese" model of ref. [Hj. The 
symmetries of these toy models allow us to find explicit solutions of the full Einstein's 
equations that again only show effects oiOipH"^), consistent with our general picture. 

The toy models will also allow us to address another problem, namely that when- 
ever inhomogeneities are present, the very definition of average cosmo logical parame- 
ters becomes ambiguous and deserves a thorough re-examination. In the literature on 
backreaction effects, the most commonly used prescription for averages is spatial av- 
erages over surfaces of constant proper time This definition is clearly not physical. 
The synchronous gauge may not even exist in a general cosmological solution, and 
even when it does it becomes singular whenever caustics in the matter flow develop. 
Moreover, it requires averaging over regions outside our past light cone. Finally, it 
gives rise to pitfalls clearly and succinctly described in 

In the toy model of Section 3 we will define several observable quantities that can 
be identified with a physical Hubble parameter. Significantly, these definitions differ 
from each other and from the homogeneous result exactly by terms 0{f'H'^). 

In Section 4 we will find arbitarily inhomogeneous cosmologies for which the Hub- 
ble flow receives no corrections at all, and briefly discuss the Swiss Cheese model, 
another inhomogeneous and anisotropic cosmology in which the local Hubble flow is 
uncorrected. 

Tellingly, the scale we flrmly associate to the large scale effect of inhomogeneities 
is O(f^), with V the typical peculiar velocity. This result is both physically sensible 
and large enough to risk to become a factor in future precision cosmology. 

Some technical material omitted in the body of the paper is collected in two 
appendices. 

2 Backreaction and Mass Renormalization 

Consider an FRW universe with metric 

ds^ = dt^ — a^{t)'jijdx'^dx^ , i,j = 1, .., 3. (1) 

In this section we will often choose a flat homogeneous space metric jij = Sij to make 
our (heuristic) equations simpler, but our results could be easily extended to other 
homogeneous metrics. 
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Linearized perturbations of this FRW universe reduce to the well-known New- 
tonian limit whenever the stress energy tensor Tmn can be decomposed into a ho- 
mogeneous piece, Tq = p, T- = 6lp, sourcing the Friedmann's equations, plus an 
inhomogeneous "dust" component to(x, t) = 6p{x,t), ~ Here an over-bar will 
denote space-averaged quantities (to simplify the discussion we will for the moment 
ignore ambiguities in the averaging proceedure, which will be addressed below). To 
first order in G6p, the metric is 

ds"^ = (1 + 2(p)dt'^ - a2(t)(l - 2(j))-f,jdx'dx^, (2) 

while obeys the Poisson equation |lj 

a-^D^D'(l) = AnG5p. (3) 

Here Di is the covariant derivative w.r.t. the metric 7ij, which is also used to raise 
and lower the indices The dominant non-linear correction to this equation has a 
simple physical meaning: gravity couples to all forms of energy, including the energy 
of the gravitational field itself. Thus, to second order in 0, equation Q becomes 

a-^D,D'(j) = 4nG6p - ^a-^D,(f)D'(j). (4) 

Though this equation is heuristic, it does capture the main effect of non-linearities 
since a more complete derivation of backreaction effects leads to a similar formula 

E]. 

Consider now the case that the background space metric is fiat 'fij = Sij, and that 
the "dust" making up 5p is composed of very compact objects of radius r distributed 
on a cubic lattice of physical size /; by compact we mean that r -C /. Then 

Sp{x) = ^ pr {ax — In) — p. (5) 

Here Pr{x) is an arbitrary positive function vanishing for |x| > r and normalized to 

j d^xpr{x) = m. (6) 

The backreaction equation (jlj) does not contain any time derivative, so, when 
studying the effect of inhomogeneities at any given time, as we will do next, we can 
set the scale factor a = 1 to avoid needless complications. 

The average background density p is a function of the cluster mass m and the 
volume of the fundamental cell of the lattice p = m/l^. 

To understand the large scale effect of the non-linear term we average eq. Q over 
a cube of side much larger than / at constant time t. The averaging procedure itself 
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introduces an ambiguity, since we could have chosen to average over a different space- 
hke surface. In this section, this ambiguity will not matter, since here we merely want 
to show how to avoid much larger, indeed divergent, unphysical effects. 
The average 5p is zero by construction, so eq. (jlj) averages to 

1. 



A0 = --V0V0, (7) 



y-+oo V Jv 



V-^oo V Jv 

The gradient and the Laplacian are the standard flat-space ones. To second order in 
G6p, the potential appearing in the right hand side of this equation is the solution 
of the Poisson eq. (jlj). 

A brief computation then shows that the average gravitational energy density is 



1 r 72 



2tt ^ 
pr \ —m 



(9) 



where pr{k) = J d^xpr{x) exp{ikx). 

Here we encounter a (pseudo) paradox: when the "clusters" are point-like, the 
distribution pr is a constant and sum in eq. © diverges. If we regulate the divergence- 
say by making Pr{x) = S/Airr^ for \x\ < r and zero otherwise-we may still end up with 
a gravitational energy density as large as the background energy Tq . This happens in 
particular if the radius of the sphere is close to its Schwarzschild radius. Physically, 
this would mean that if all matter in a universe at critical density was distributed 
into black holes, corrections to the linear approximation would be 0(1) no matter how 
small the lattice step b\. Even when r 3> Gm this contribution is suspicious because 
it depends strongly on the size of the cluster. 

An estimate of the r dependence of the sum is given in the Appendix, here we 
only quote the result, and we give a simple argument to justify it: 



SttG 



V0V0 = -Cm' 



rP l^ V/5 



(10) 



where Ci, C2 are dimensionless numbers of order unity which are determined by 
the density profile of compact objects pr and by the shape of the lattice. The first 
term in the expansion (fTIH) comes from the classical gravitational energy of a body 
of mass m and size r. Other terms, describing corrections, must be there because for 
a body of uniform density filling the entire lattice cell t^ vanishes by definition. The 
divergence at small r term scales with / as pressureless mass density, while the finite 
term, independent of r, scales as ultra-relativistic energy density. 
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At this point, it is obvious that the divergent term is unphysical and it must be 
canceled by redefining the "bare" mass of the compact objects to first order in G by 

m = ruphysicai + Gi—^. (11) 

This classical renormahzation is not an option because the gravitational field is actu- 
ally determined by the physical mass. 

- Y / d^Xpr{Q)^^pr{x) =m-G J ^^^^ Prik) ^ = niphysical- (12) 



Notice that corrections due to inhomogeneities do exist. They are due to the finite 
term in eq. (jlUp . which cannot be eliminated by a renormahzation of the mass, since 
it scales with / as radiation. These corrections change the Friedmann equation by 
terms 0{H^P). If / is interpreted as the non- linear length scale of the actual universe, 
then HH'^ ~ f , where v is the typical peculiar velocity. So, non-linear corrections to 
the Hubble parameter due to gravitational energy are of the same order as relativistic 
corrections due to pecuhar velocities. We have actually neglected the latter as well 
as ambiguities in the definition of the Hubble parameter. Evidently, to make further 
progress we need to be more systematic, even if at the price of studying a vastly sim- 
plified model of inhomogeneities. This is what we will do next, starting by considering 
a distribution of pressureless matter that only breaks translational invariance in one 
direction, while preserving translations and rotation in two orthogonal directions. 



3 Dust- Wall Universe 

Our goal is to construct an exact solution of the Einstein's equations describing an 
expanding inhomogeneous universe. A simple possibility is planar symmetry - the 
metric depending on t and x but not on y and z, and isotropic in the yz plane. Here 
we use Taub's jS] explicit expressions for the metric of vacuum plane-symmetric space- 
times to construct a universe of equidistant plane-parallel dust walls. Our treatment 
of the walls is similar to the single wall case fOJ • 
The plane-symmetric metric can be written as 

ds^ = e2"(rft2 -dx^)- e^^'idy^ + dz^), (13) 

where u = u{t, x) and v = v{t,x). The walls are located at x = ±b, ± 36, ... so that 
the metric coefficients u and v are 26-periodic functions of x. It also follows from the 
symmetry of the problem that u and v are even functions of x. 
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The stress-energy tensor of a thin dust wall of proper surface density a{t) located 
at X = 6 is 

Ti;{t, x) = (T(t)e-"(*'^)(5(x - b) diag(l, 0, 0, 0). (14) 

For simplicity we consider walls with zero pressure, but it is easy to find the analogous 
solutions for an array of walls with general equation of state. Calculating the Einstein 
tensor for the metric (fT^. equating it to the stress-energy tensor of the dust walls, 
and using the symmetry of the metric, we get the following jump conditions at the 
X = b wall: the (jQ-component gives 

a = 4e-"9^w, (15) 

the Gl and G\ components are non-singular at the wall, and the = components 
give 

d,{u + v) = 0, (16) 

where the derivatives are calculated for x = b — 0. 

Between the walls the metric is given by the Taub's expressions: 

^2n^ /V ^ ^2v^f^g^ f = f^t + x), g = g{t~x), (17) 
Vf + 9 

where / and g are arbitrary functions of one variable, and /' and g' are their deriva- 
tives. Since u and v are even functions of x at all time, we must take g{^) = /(O- 
Then (jlfij) gives an equation for /: 

fit + b) fit - b) 1 fit + b)- fit -b) _ 

f>(t + b) f'it-b) ^2 f{t + b) + f{t-b) • ^ ^ 

The remaining jump condition (fTH|l gives the proper surface density as a function of 
time in terms of /. 

At large times t^b there are many walls within the horizon and we recover the 
homogeneous matter-dominated universe. To lowest order in 6, finite differences in 
equation ([T^ can be replaced by differentials, and we get 

(f)Vif^O. (18) 

with the solution f{t) oc t^. To lowest order in b this gives the metric e^" oc e^^ oc t^. 
This describes homogeneous matter-dominated universe with conformal time t. 

As expected, there are finite b corrections. In the next to leading order, we get 
f{t) oc + AOb^t'^. This gives the metric coefficients 

oc + f{20b^ - 6x2), ^2v ^ ^4 ^ ^2^4Q^2 ^ Q^2y (20) 
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We see that the metric is close to that of a uniform matter dominated universe, with 
h"^ /t^ corrections. In physical units the corrections are ~ HH"^, where I is the physical 
distance between the walls and H is the Hubble constant. 

Now wc would like to see how the averaged expansion rate of the universe is 
changed due to inhomogeneities. To this end, we define parallel and perpendicular 
scale factors on the wall: 

a\\{t) e<*^''\ ax(t) oc e"(*'^). (21) 

In the next to leading order we get 

a\\ (t) oc + 236^, a^{t) (xt^ + 76^. (22) 

Further we define parallel and perpendicular Hubble constants 

d\oga\\ dloga± 

where dr oc a±dt is the proper time on the wall. Both definitions of the Hubble 
constants are physical. H\\ can be measured from the proper surface energy density 
as a function of proper time. H± can be measured by light propagation time be- 
tween the walls. Wc then get corrections to the Priedmann equation for parallel and 
perpendicular Hubble constants 

H^al oc 1 + ^HH\ Hlal a 1 - ^HH', (24) 

where I 2be" is the physical distance between the walls and H Ki H\\ H± 2/ (3t). 



4 Uncorrected Hubble Flow 

In this section we will consider a broad class of examples of anisotropic cosmologies 
for which wc can find exact solutions, this time sourced by arrangements of parallel 
relativistic cosmic strings rather than parallel walls. The dynamics of such objects- 
so long as they remain parallel-is identical to that of massive point particles in 2+1 
dimensional gravity, and so for the rest of the section we will use that language. 
We will show that starting with any isotropic and homogenous solution of Einstein's 
equations in 2+1 dimensions one can add point particles-or indeed an arbitrary dis- 
tribution of dust-and (due to the special properties of gravity in 2+1 dimensions) the 
presence of the point particles and/or dust does not modify the Hubble expansion at 
all. 



7 



We begin with the metric 

^ _ a{tyds2^ (25) 

where as usual ds^ is the metric of a homogeneous and isotropic two-dimensional flat, 
spherical, or hyperbolic space. Prom this ansatz the 2+1 Priedmann equation follows: 

H^^p-k/a^, (26) 

where as usual k — {—1,0,+!}, H = a/a and p = pQ{ao/ ay^^'^'^^ for the case of a 
perfect fluid with equation of state wP' 

One peculiar feature special to the case of matter domination [w — 0) is that 
a{t) — t, with zero acceleration. This indicates that uniformly expanding dust in 2+1 
dimensions does not exert a force on itself, a result which might have been expected 
since there is no force between static particles in 2+1 dimensions. Note also that, 
since in this case the matter density scales the same way as curvature, the scale factor 
is the same for all three types of universes. 

4.1 Point Particles 

In order to consider an anisotropic cosmology we begin by reviewing the physics of 
a point particle with mass m G (0,27r) in a 2+1 dimensional space'^. As is well- 
known the metric is simply flat space with a conical defect along the world-line of the 
particle: dsl = dt^ — dr"^ — r'^dO'^, where the range of 9 is from to 27r — m. Perhaps 
the simplest way to see this is to recall that in 2+1 dimensions Einstein's equation 
in vacuum R/^i, — implies R^vXa = 0, so that the space must be flat away from 
the point particle. The conical defect induces a delta function in the curvature of 
strength 2m, from which it follows that the metric satisfles Einstein's equations with 
a point-like source. Since any space is locally flat, all of the solutions we consider 
here are of this form very close to the singularity. Note that the mass of the particle 
should be smaller than 27r since the circumference of the circle goes to zero in that 
limit. 

Let us consider a collection of point particles with masses at locations Zi in 
flat space. Taking the metric ansatz 

refers to the 2+1 dimensional equation of state. Lifting w = here to 3+1 gives an anisotropic 
but homogeneous 3+1 dimensional universe filled with a uniform "gas" of parallel cosmic strings. 
•^We will work in units where SttG^ = 1. 
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Einstein's equations reduce to 



N 

d,d,(l>=-Y.m,5\z-z,). (28) 

i=l 

So long as the sum of the masses satisfies the constraint M = Y^iLi < 27r the 
solution is 

ef = Ill,\z-z,\-'^. (29) 

It is easy to see where the constraint on the total mass comes from. Consider a large 
disk S containing all the masses. Recall that the Euler character x = 1 for ^ disk, 
that the space is fiat away from the point particles, and that a point particle of mass 
m gives rise to a delta function in the scalar curvature R of strength 2m. Then, since 
the integrated extrinsic curvature of a circle is the derivative of its circumference with 
respect to the radius, drC, the Gauss-Bonnet theorem. 



R + 2 K = 47rx, (30) 
s JdS 

implies = 27r — 9^0. If the total mass M exceeds 27r the space cannot be globally 
fiat, since the circumference of the disk must shrink as the radius increases. If we 
insist that the curvature be zero, the solution (j^^ now implies the presence of a 
conical singularity at z = cxd equal to Air — Y,f=i ^i- The metric then describes a 
compact space with spherical topology, with + 1 conical singularities whose total 
mass adds up to An. 

From the Gauss-Bonnet theorem it follows that in order to accommodate point 
particles with total mass M > An some curvature is needed. We wish to consider 
the case of a large or infinite number of (positive) masses, and so the background 
2-dimensional curvature must be negative. In fact the expansion of the universe 
automatically acts as a source for the negative spatial curvature. 

Consider a metric of the form 

ds^ = dt^ - a{tfe'l'^'~'^dzdz. (31) 

If we now require that the energy-momentum tensor consists of point particles of 
mass rrii at fixed spatial locations Zi plus a homogeneous background matter density 
Phom, Einstein's equations reduce to 

k ^ 

dzdgcp = --e'f' -^niiS'^iz - Zi), 



k 

— = Phom 5", (32) 

a) 
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where k determines the (constant) spatial curvature away from the singularities. This 
is the main result of this section. As is clear from the second equation the evolution 
of the scale factor is not disturbed by the defects as long as the first equation can be 
satisfied. This in turn depends only on the masses rrij (see appendix). The conclusions 
above hold even in the presence of general homogenous matter components. In this 
case the pressure of the homogenous background is given by, 

P = --. (33) 
a 

For the case A; = the first equation reduces to the Poisson equation, whose 
solution with point particle masses we discussed above. When k = ±1 eq. 
is the Liouville equation with positive and negative curvature respectively. As we 
discuss in the appendix, solutions to this equation always exist, unless topological 
obstructions arise. From the Liouville equation it also follows that, as required by 
the Gauss-Bonnet theorem, the local density of particles cannot be too large to avoid 
"overdosing" the universe. 

One case of interest is a square grid with a point mass at each vertex and with the 
background matter density phom = 0. Since the total mass in the universe is greater 
than Att the background curvature must be negative. The mathematical problem 
we need to solve reduces to finding the metric with constant negative curvature and 
a single conical singularity on a square torus. Although in this case the Liouville 
equation cannot be solved in closed form, the solution is known to exist and to be 
unique. Physically, while the metric has locally negative spatial curvature away from 
the particles, it has precisely the correct density of positively curved defects to cancel 
the negative curvature when averaged over a large region (or over one tile of the grid). 
The solution approximates flat FRW at large scales (and it is for this reason that it 
is possible to arrange the masses in a square grid even though the metric locally is 
negatively curved). 

In general there is a very large continuous infinity of solutions of this type (cor- 
responding to the freedom of changing the positions and masses of the particles). In 
fact (as is clear since we can take the limite of very small strength 5-functions closely 
spaced) if we replace the term Y^iLi'^^i^'^i^ ~ ^i) equation (jH^ with a function 
Pd\ist{z,z) we can obtain solutions with arbitrary distributions of dust. 

It is remarkable that the function a{t) in the metric ()25|) remains unchanged 
for any arrangement of point masses or distribution of dust, and that therefore the 
local Hubble expansion law is entirely uncorrected by the inhomogeneities. This 
is a physical statement confirmable by local experiments; for example an observer 
in this universe could measure the rate of change in the local matter density as 
a function of his proper time and would observe it to be determined only by the 
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average Hubble constant (averages here can be taken over surfaces of constant t, for 
which the expansion is uniform), independent of his location and the locations of the 
anisotropies. 

Let us mention that, even though the Hubble expansion is not modified, many 
experiments will be sensitive to the positions of the masses, for example test particles 
fired out along some trajectory will scatter off the point masses. For the same reason 
the observed brightness of a distribution of standard candles of fixed red-shift will 
vary. 

4.2 Swiss Cheese 

There exists in the literature another exact solution describing an inhomogeneous 
universe which gives no corrections to the Hubble parameter-the "Swiss Cheese" 
model of Schucking |n|. This is a universe in which spherical regions of homogeneous 
background matter are excised and replaced with vacuum, but with a spherical con- 
centration of mass at the center of the void. If the mass at the center is chosen 
appropriately and the voids are non-intersecting the metric outside the voids is com- 
pletely unaffected by their existence and the Hubble flow as measured by an (outside) 
observer is unchanged. Of course as before there are many experiments which would 
be sensitive to the inhomogeneities. 

5 Conclusion 

Both generic (but heuristic) as well as exact (but highly symmetric) models show 
that the effect of matter inhomogeneities on the large scale expansion of the universe 
is small, bounded from above by v"^ ~ 10~^. The numerical coefficient and sign of 
the effect depend on the model and the definition chosen for Hubble. What we have 
shown can also be interpreted as follows: there always exist gauges where the metric 
differs from the FRW metric only at order ~ f^. 
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Appendix A: Estimating the Gravitational Self-Energ 
Divergence 

Here we estimate the behavior of the sum in eq. Q- 

The solution of the Poisson eq. (j3]) is easily written in Fourier transform as 



1,2 ^ /3 

mGZ3,m^0 



27r 

T 



(A.i; 



For I —>■ oo, eq. Q can be approximated by the integral 



hm — 



E 



27r 



-m 



d^k 2 



TC 



(27r)3 F 



Pr{k) 



(A.2) 



Since Pr{k) is analytic and square summable, we can also estimate the sum at finite 
I as 

1 ^ f /2n 



/3 E 27rr772 



9 o 

m /r 
" ~^ [l 



(A.3) 



By dimensional analysis, f{x) is a dimensionless function, independent of m,r,l, and 
smooth in a neighborhood of x = 0. 

Therefore we can see that the energy density diverges as 1/r for r — > 0, and we 
also recover the finite term of eq. (fTUj) . 



Appendix B: More on the Liouville Equation 

As we have shown in section IH the problem of point particles (with zero peculiar 
velocities) in an homogenous FRW universe in 2+1 dimensions reduces to finding 
solutions to the Liouville equation, 

k ^ 
^ i=i 

This equation describes a two dimensional surface of constant curvature k with conical 
singularities m^. The general solution of the Liouville equation is 

= (B.2) 
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where w{z) is a holomorphic function that must be chosen appropriately in order to 
reproduce the correct singularities at 



~ -log\z — Zi\, as z ^ Zi. (B.3) 
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In principle w{z) can be determined by using the technology of the Fuchsian equa- 
tions (see for example [13) 5 even though explicit solutions can only be found in 
special cases. The solutions of the Liouville equation are in general determined by 
the singularities and by the topology of the space. Since the locations Zi of the singu- 
larities in eq. ()B.1|) are arbitrary (up to conformal transformations) there is actually 
a large moduli space of solutions. Note that even though the Zj's are free, the posi- 
tion in physical space are constrained since on scales smaller than the curvature the 
total deficit angle cannot exceed 27i. This is automatically guaranteed by solving the 
Liouville equation. 

Let us now consider the case where the space is compact. By integrating both 
sides of eq. ()B.1|) one finds, 

V=Mi^4^S£^, (B.4) 

k 

where is the volume and (7 is the genus of the surface. This formula follows from the 
fact that the left hand side of eq. ()B.1|) is proportional to ^/gR and the Gauss-Bonnet 
formula for compact surfaces without boundaries. 



v^i? = 2 - 2g. (B.5) 



1 

4:71 Js 

Since the volume must be positive, eq. ()B.4|1 implies restrictions on the allowed 
singularities depending on the sign of the curvature. 



k = l 

This situation can arise for example in the presence of background dust or positive 
cosmo logical constant. Since we only allow positive masses the volume formula ()B.4|1 
here determines that the only possible topology is the spherical one and that the 
sum of the deficit angles must be less than 47i. The Liouville equation can be solved 
explicitly in the case of three singularities. In general a solution exists provided that 
a mild constraint on the deficit angles is satisfied, namely that the largest conical 
defect is less than the sum of the remaining ones jTHj . 

■*This asymptotic behavior holds for rrii < 2tt. When k is negative singularities with nii = 2tt (so 
called parabolic) are also acceptable in which case ~ — 2 log \z — Zi\ — 2 log | log |z — Zi| |. 
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k = -l 



This is the case relevant for particles in vacuum when the sum of the deficit angles is 
greater than An. An important theorem dating back to Poincare and Picard estab- 
lishes that when the curvature is negative a solution of the Liouville equation with 
prescribed singularities always exists and is unique, once eq. ()B.4|) is satisfied. If the 
sum of the deficit angles is greater than An any topology is allowed and in particular 
we can choose the spherical topology. Note that in hyperbolic space a deficit angle 
increases the volume of the space. A case relevant for our purposes is the torus with 
one conical singularity. This problem is equivalent to that of a grid of equal masses. 
From eq. ()B.4|) we derive that the volume of each tile equals the deficit angle. This 
ensures that the average curvature is zero as expected on physical grounds. 
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